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SAR 图 像 配 准 的 二 元 Rayleigh 分 布 互 信息 方法 * 
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摘 X 合成 孔径 雷达 (Synthetic Aperture Radar, WIR SAR) 图 像 中 乘 性 斑点 噪声 的 存在 ， 使 得 传统 
的 互信 息 配 准 方法 中 的 插值 假象 更 为 严重 . 本 文 把 单 视 SAR. 幅 值 图 像 同 质 区 域 服 从 Rayleigh 分 
布 这 一 先 验 信 息 引 入 互信 息 配 准 方法 中 ， 对 同 源 单 视 SAR 幅 值 图 像 进 行 配 准 . 以 合成 图 像 和 真 
实 图 像 的 试验 结果 和 分 析 证 实 了 方法 的 可 行 性 ， 说 明 新 方法 优 于 传统 直方 图 估计 联合 分 布 的 互信 


息 方 法 . 
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1 引言 


基于 互信 息 的 配 准 方法 具有 全 自动、 精度 高 、 无 需 提 取 图 像 特征 等 优点 ， 在 医学 图 像 配 
准 丰 和 遥感 图 像 配 准 2-4 中 得 到 了 广泛 的 应 用 . 但 对 于 SAR 图 像 ， 由 于 其 成 像 机 理 决定 的 乘 
性 斑点 噪声 的 存在 ， 使 得 互信 息 配 准 方法 中 的 插值 假象 问题 更 为 突出 ， 从 而 限制 了 该 方法 
E SAR 图像 配 准 中 的 应 用 . 文献 [5,6] 都 是 先 对 SAR .图像 进行 滤波 预 处 理 ， 再 利用 互信 息 实 现 
了 SAR 图像 的 配 准 . 

传统 的 互信 息 配 准 方法 不 考虑 图 像 的 任何 先 验 信息 ， 只 用 统计 直方 图 估计 联合 分 布 来 计算 
互信 息 . 而 文献 [7] HRH T SAR 能量 图 像 服从 而 Gamma 分 布 这 一 先 验 信 息 ， 根 据 单 幅 图 像 
的 分 布 类 型 在 两 幅 图 像 之 间 建 立 二 元 Gamma 分 布 模型 ， 然 后 利用 图 像 数 据 估 计 模 型 参数 进而 
计算 互信 息 ， 较 传统 的 互信 息 方法 更 准确 . 但 SAR 图 像 常 以 幅 值 的 形式 给 出 ， 文 献 [8-10] 说 明 
了 单 视 SAR. 幅 值 图 像 的 同 质 区 域 服从 Rayleigh 分 布 . 

本 文 在 上 述 文 献 的 基础 上 对 两 幅 单 视 SAR. 幅 值 图 像 建立 二 元 Rayleigh 分 布 模型 ， 由 二 
元 Rayleigh 分 布 的 密度 函数 得 到 了 二 元 Rayleigh 分布 的 互信 息 的 形式 和 计算 方法 : 给 出 了 
两 幅 单 视 SAR 幅 值 图 像 配 准 的 二 元 Rayleigh 互信 息 配 准 算法 . 利用 该 算法 对 合成 图 像 和 真 
实 SAR 图 像 进 行 试 验 ， 并 与 传统 的 互信 息 方法 作 了 比较 和 分 析 ， 结 果 表 明 该 方法 可 用 于 同 源 
单 视 SAR 幅 值 图 像 的 配 准 ， 较 传统 互信 息 方 法 更 有 效 . 


2 二 元 Rayleigh 分 布 及 其 互信 息 


本 节 介 绍 二 元 Rayleigh 分 布 及 其 互信 息 ，2.1 节 给 出 了 二 元 Rayleigh 分 布 及 其 参数 估 
计 ，2.2 节 给 出 了 二 元 Rayleigh 分 布 互 信息 的 推导 过 程 ，2.3 节 给 出 了 二 元 Rayleigh 分 布 的 
模拟 数据 产生 方法 . 
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2.1 二 元 Rayleigh 分 布 及 其 参数 估计 
二 元 Rayleigh 分 布 (Bivariate Rayleigh Distribution, BRD) 是 Nakagami-m 分 布 在 m= 二 1 时 
的 特殊 情形 0， 密度 函数 为 


47129 ( 1 (2 pa ) 2 /Dx122 
ajo mm as(- (Huy E), a 
Ta xs (2,22) a3a2(1 — p) SAC gesta: i az /) \(1 — p) /a1a2 Y) 


其 中 
Tı > 0, T2 > 0, Ql = E(X1), G2 = E(X7), 
Cov(X?, X2) 2 1 z Ņ\ 2k 
EE A IN a A N E 
í Var(X?) VVar(X2) oi) >, kIT (k +1) ( 2 ) 
E EN — RJE Bessel FRA. 
二 元 Rayleigh 分 布 的 边缘 分 布 为 
2 
fx (0) - exp ( - -), y »O.. wq. (2) 
iX EAE JL Rayleigh 分 布 的 密度 图 数 ， 也 是 (1) 式 之 所 以 为 二 元 Rayleigh 分布 的 原因 之 一 ， 也 
是 可 以 用 它 对 同 源 的 两 幅 SAR 幅 值 图 像 进行 建 模 的 基础 . 
二 元 Rayleigh 分 布 的 各 阶 矩 为 





m+n mMm mn z 1 
E(X"Y")—(1— p) t F afa? D TR+t1l+t E )r(e+1 +2), (3) 


HP m, n 为 非 负 整数 .从 而 得 到 期 望 和 方差 








V Tas 4 一 
Exo i: Var(X;)— ———-a;, i— 1,2, 


以 及 (1) 式 中 各 参数 的 意义 . 
给 定 样本 (21, 723), (01,23). (mp, 2) 后， 可 用 和 矩 估计 来 估计 各 参数 


20 l« | : > 
ái = m 1 — 1,2, pom - : * 
= È (etr E P-a) 


2.2 ZIT Rayleigh 分 布 的 互信 息 
用 密度 函数 表示 的 互信 息 的 形式 为 


T. py(z, y) 
B)= | f funr) m Tee dad, (9 
其 中 feas (e, y) 为 随机 变量 A 和 8B 联合 密度 函数 函数 , 为 联合 密度 函数 函数 的 文 撑 集 , falx) 
和 fp(zx) 分 别 为 A4 和 8B 的 边缘 分 布 . 由 (4) 互信 息 还 可 以 表示 为 


f(A4,B) (A, B) ) 


MI(A, B) = E(In IESU 


(5) 
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其 中 E(-) 是 对 联合 分 布 fia,B) (2, y) 取 期 望 . 

先 说 明 几 个 记号 ， 由 (1) 确定 的 二 元 Rayleigh 分 布 可 由 三 个 参数 a1, a3 和 p 来 确定 ， 因 此 记 
之 为 BRD(ai,a2,p)， 对 应 的 互信 息 记 为 MI(a1, az,p)， 由 于 


n {CX XX) _ mn (i) tee (n (PEST )) - (9L ) 





fx (X) fx (X2) 1 — p)J/a1a2 l~p\a a2 
于 是 由 (5) 式 可 得 
加 fo xy OX, X2) 
MI(a1,a2,p) = E(In P 39) 
c D ür 2) EUER E(In (olg m en) 


HF EC) ZEE 218 BRD(a1,a2, p) RWE, id 


1X2 
A - (mn (c2) 


这 个 积分 不 能 给 出 显示 表示 ， 但 是 进一步 有 








A = FE(m Ra 
z f |a TEL -— retta) (CX) In (r(A) azaz 
= Ë(m e 


Hp EC) 是 对 BRD(1, 1,p) RHE. 最 后 得 到 的 积分 还 是 不 能 得 到 显示 的 表达 式 ， 但 它 已 经 是 
参数 p 的 函数 ， 即 只 与 一 个 参数 有 关 ， 从 而 MI(ai, a2, p) 可 记 为 MI(p). 在 实际 应 用 中 给 定 参 
数 p， 用 数值 方法 给 出 数值 结果 . 通过 数值 计算 ， 可 以 得 到 MI(p) 和 pp 的 关系 如 图 1 所 示 . 


MI(p) 


图 1: 二 元 Rayleigh 分 布 的 互信 息 MI(p) 和 其 参数 p 的 关系 
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2.3 二 元 Rayleigh 分 布 的 模拟 数据 产生 方法 
为 了 分 析 二 元 Rayleigh 分 布 的 互信 息 在 图 像 配 准 中 的 可 行 性 ， 首 先 作 合成 图 像 的 配 准 试 
验 ， 从 而 要 产生 服从 二 元 Rayleigh 分 布 的 模拟 数据 . 由 (1) 可 得 的 边缘 分 布 的 分 布 函 数 为 


2 
Fx, (xı) = 1 — exp =R >0, (6) 
^ Fx, (21) = W 4 ul 为 [0， 1] 上 的 均匀 分 布 ， 则 


ĉı V —ajln(1-— ui) 


为 服从 (6) 式 的 采样 值 ， 同 样 由 (1) 可 以 得 到 条 件 分 布 


272 pazĉ? + air? 2VDZ172 
X= = —— (- Za) (————), £ > 0, 
Papa (1X = £1) aia2(1 — p) m aja5(1— p) ^ ^N /a1àz(1 — p) i 


从 而 


2 2 
Pam ll n=) = 1-A (mfr yaa ) 


HEF Qla, b) X Marcum-Q AG, EXN 
OO rx? +a? 
dps i kaniz 
Q(a, b) f zexp( : ) o(az)da 


S Fyx (22X1 = ĉi) = us; F ua A [0,1 上 的 均匀 分 布 ， 则 解 该 方程 得 的 $2 为 服从 的 分 
HL. FÆ EZ Rayleigh 分 布 的 步骤 如 下 : 
1y up 
ĉi = V/—ajlni(1— ui), 3X $i y-a lnu, 


其 中 ui 为 [0,1] 均匀 分 布 的 采样 值 ; 
2) IHR ĉo MAJ fz Q(a,t)-1-— us EK Q(a,t) = uo, H uz X [0,1] 均匀 分 布 的 采样 值 ， 且 


m DT RT Vas 2i 


解 得 的 (21,22) 就 是 来 目 BRD(al, a2, p) 的 一 个 模拟 数据 . 





3 SAR 图 像 配 准 的 二 元 Rayleigh 分 布 的 互信 息 方 法 步骤 


给 定 竺 配 准 的 两 幅 SAR 图像， 建立 二 元 Rayleigh 分 布 模型 如 下 : 给 定 一 个 变换 ， 把 参考 图 
ee ban wi AER —— 


NM AD CENT MICA B. 

TEBOHEZ Bj. Bu E RHR RIAA N AE Xe f. ROS e Be RE. SLT Be DR uE EN] 
服从 的 分 布 类 型 是 Rayleigh 分 布 ， 在 配 准 的 过 程 中 ， 所 选 定 的 区 域 如 果 不 是 同 质 区 域 ， 即 使 分 
布 类 型 是 Rayleigh 分 布 ， 但 参数 不 同 ， 佑 计 参 数 和 计算 互信 息 时 也 是 有 问题 的 .因此 在 用 该 方 
法 配 准 之 前 ， 先 用 Kolmogrov 检验 选 定 区 域 . 





Kolmogrov 检验 不 仅 需要 知道 分 布 的 类 型 ， 还 需要 知道 分 布 的 参数 . 因此 先 用 Rayleigh 分 
布 的 窍 估 计 方 法 估计 出 参数 ， 再 用 Kolmogrov 检验 方法 检验 样本 是 否 为 来 自 该 参数 的 Rayleigh 
分 布 . 

”通过 Kolmogrov 检验 满足 分 布 要 求 的 区 域 给 予 保留 ， 不 满足 则 舍弃 ， 重 新 选取 区 域 进行 检 
验 ， 要 是 选取 多 次 都 被 拒绝 ， 可 能 是 选取 的 区 域 非 同 质 ， 也 可 能 图 像 本 身 就 不 符合 要 求 . 对 于 
前 一 种 情况 可 以 继续 重新 选取 区 域 或 人 工 选 取 或 先 粗 略 分 割 图 像 ， 阁 是 后 一 种 情况 停止 用 该 方 
法 配 准 . 

单 视 SAR 幅 值 图 像 的 二 元 Rayleigh 分 布 互信 息 配 准 方法 步骤 如 下 : 

1) 记 参 考 图 像 和 竺 配 准 图 像 为 和 sf Men HERE- ARE To, A TEREFE 
快 和 优化 时 不 陷入 局 部 极 大 值 ， 最 好 给 定 的 初始 变换 能 粗略 配 准 图 像 ， 例 如 可 用 人 工 选 取 控 制 
点 的 配 准 方法 确定 初始 变换 ，Isen 变换 后 得 到 的 图 像 记 为 hrmt(To); 

2) 用 Kolmogrov 检 验 在 Ler 和 mr(7T6) 的 重合 部 分 选取 用 才 估 计 参 数 和 计算 互信 息 的 区 
域 ， 选 取 区 域 可 以 随机 选取 也 可 以 人 工 选取 ， 为 了 较 快 选 定 满 是 要 求 的 区 域 ， 人 工 选 取 同 质 区 
域 较 好 ， 也 可 用 分 割 的 方法 选取 ， 当 然 这 里 对 分 割 的 要 求 不 高 只 需 把 同 质 区 域 分 割 出 来 即 可 ， 
所 选区 域 可 是 一 个 也 可 为 多 个 ; 

3) H Powel 的 方 回 集 优 化 算法 由 初始 变换 1 最 大 化 互信 息 得 到 最 优 变换 思 nt， 其 中 计算 
互信 息 的 过 程 如 下 ， 在 无 sf 中 保留 由 2) 步 选 定 的 区 域 ， 对 于 每 给 定 一 个 变换 荆 和 得 到 的 变换 图 
像 frmf(T)， 用 Lec PETA EKERI hem (D 中 的 重合 部 分 估计 二 元 Rayleigh 分 布 的 参数 ， 
进而 计算 互信 息 MI(T)， 需 要 注意 的 是 如 果 选 定 的 区 域 为 多 个 ， 这 多 个 区 域 并 不 一 定 是 来 自 同 
一 个 二 元 Rayleigh 分 布 的 样本 ， 因 此 对 每 一 个 区 域 都 计算 ~ 个 互信 息 值 ， 最 后 的 MI(T) 为 各 区 
域 互 信息 值 之 和 . 

得 到 最 优 变 换 工 wt， 便 完 成 了 两 幅 图 像 的 配 准 .本 文 用 的 是 仿 射 变换 ， 根 据 不 同 图 像 的 要 
求 ， 也 可 选 其 它 的 变换 . 值得 注意 的 是 在 变换 过 程 中 所 用 的 插值 方法 也 很 重要 ， 本 文选 用 双 线 
性 插值 . 





4 ”试验 结果 与 分 析 


我 们 分 别 用 合成 图 像 和 真实 SAR 图 像 做 了 配 准 试验 ， 并 对 结果 进行 了 比较 和 分 析 . 

4.1 合成 图 像 实 试验 分 析 

图 2 是 一 组 服从 二 元 Rayleigh 分 布 的 合成 图 像 ， 左 边 两 幅 图 像 大 小 完全 相同 ， 区 域 边 界 划 
分 完全 一 致 ， 同 一 位 置 的 两 个 像素 值 是 一 个 米 自 二 元 Rayleigh 分 布 的 模拟 数据 . 两 幅 图 像 各 有 
四 个 不 同 的 同 质 区 域 ， 四 个 区 域 参 数 分 别 是 (25, 27,0.9), (45, 50, 0.8), (56, 62, 0.9) 和 (70, 71,0.9). 
右边 两 幅 图 像 是 从 左边 两 幅 图 像 中 截取 的 部 分 图 像 ， 但 是 截取 的 区 域 不 同 ， 从 而 这 两 幅 图 像 有 
重 营 部 分 但 又 不 完全 相同 ， 可 用 于 图 像 配 准 试验 分 析 ， 需 要 的 变换 只 是 平移 . 








图 2: ”两 幅 服从 二 元 Rayleigh 分 布 的 合成 图 像 
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把 图 2 第 三 幅 图 像 作为 参考 图 像 ， 第 四 幅 图 像 作 为 待 配 准 图 像 ， 准 确 配 准 结果 为 待 配 准 图 
像 向 左 和 向 下 各 平移 56 个 像素 ， 如 图 3 左 图 所 示 . 图 3 中 图 和 右 图 分 别 为 直方 图 估计 的 互信 息 
和 基于 二 元 Rayleigh 分 布 的 互信 息 两 个 测度 在 最 优 配 准 结果 两 个 平移 方向 附近 的 值 . 从 这 两 个 
图 中 可 以 明显 看 出 直方 图 估计 的 互信 息 振 荡 很 大 ， 不 能 用 于 图 2 所 示 的 两 幅 图 像 的 配 准 ， 而 基 
T Rayleigh 分 布 的 互信 息 测度 可 以 衡量 配 准 的 精确 程度 . 
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图 3: 最 优 配 准 拼接 结果 和 最 优 配 准 附近 的 两 种 互信 息 值 


4.2 SAR 图 像 试 验 
图 4 为 两 幅 同 一 区 域 不 同时 刻 的 SAR 图像， 图 像 大 小 均 为 512 x 512. 试验 中 确定 左 图 为 参 
考 图 像 ， 右 为 待 配 准 图 像 ， 选 用 仿 射 变换 . 





图 4: 同一 区 域 不 同时 刻 的 两 幅 SAR 图 像 


图 5 为 图 4 的 配 准 拼接 结果 ， 左 图 为 人 工 选取 4 个 控制 点 得 到 的 配 准 结果 ， 中 图 和 右 图 分 别 
是 用 传统 的 直方 图 估计 的 互信 息 方法 和 基于 二 元 Rayleigh 分布 的 互信 息 方法 由 初始 变换 出 发 得 
到 的 配 准 结果 . 配 准 结果 的 显示 方式 是 对 于 重合 区 域 各 显示 参考 图 像 和 变换 图 像 的 50%， 对 空 
白 区 域 即 参考 图 像 和 变换 图 像 都 不 能 覆 盖 的 像素 位 置 显示 为 黑色 . 从 图 5 中 可 以 看 出 ， 左 中 有 
许多 明显 的 线条 都 没有 对 齐 ， 中 图 的 结果 较 左 图 要 好 ， 但 是 有 些 线条 和 亮点 都 显得 非常 模糊 即 
配 准 结果 不 够 精确 ， 但 基于 二 元 Rayleigh 分 布 的 互信 息 方 法 配 准 结 末 右 图 中 重合 部 分 线条 清 
晰 ， 参 考 图 像 和 配 准 图 像 中 的 强 散 射 亮点 也 对 齐 的 很 好 ， 总 体 效 果 令 人 满意 . 

图 6 为 两 幅 同 一 区 域 不 同 尺 度 的 SAR 图 像 ， 这 两 幅 图 像 尺 上 度 相 计较 大 ， 左 图 大 小 为 300 x 
300. ARX 150 x 150. 以 左 图 为 参考 图 像 ， 右 图 为 待 配 准 图 像 ， 仍 选用 仿 射 变换 .图 7 为 配 
准 结果 ， 配 准 结果 的 显示 方式 同 图 5， 从 中 可 以 看 出 ， 左 图 中 细 长 的 特征 由 于 没有 对 齐 显示 成 
了 两 条 ， 重 合 部 分 到 非 重合 部 分 的 过 渡 非 常 不 自然 ， 中 图 是 由 左 图 互信 息 优化 得 到 的 ， 虽 然 把 
左 图 中 的 一 些 没 有 对 齐 的 特征 对 齐 了 ， 但 又 多 了 一 些 其 它 特征 没有 对 齐 ， 右 图 结果 则 明显 优 于 
左 图 和 中 图 . 
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图 5: 图 4 配 准 拼接 结果 





图 6: 同一 区 域 不 同 尺度 的 两 幅 SAR .图像 





图 7: 图 6 配 准 拼接 结果 


5 ”结论 


本 文 由 二 元 Rayleigh 分 布 的 密度 函数 推导 了 二 元 Rayleigh 分 布 的 互信 息 ， 对 服从 Rayleigh 
分 布 模型 的 符 配 准 的 两 幅 同 源 单 视 SAR 幅 值 图 像 建立 了 二 元 Rayleigh 分 布 模型 ， 并 提出 了 这 
样 两 幅 图 像 的 二 元 Rayleigh 分 布 的 互信 息 的 配 准 方法 . 合成 图 像 和 真实 图 像 的 试验 结果 证 实 了 
方法 的 可 行 性 ， 并 较 传 统 的 互信 息 方 法 有 效 . 

文中 对 二 元 Rayleigh 分 布 的 参数 估计 采用 较为 简单 的 矩 估 计 ， 这 个 估计 方法 可 能 不 够 准确 
从 而 会 影响 配 准 结果 ; 55h, —7u Rayleigh 分 布 的 互信 息 的 最 后 形式 虽然 只 和 一 个 参数 有 关 ， 
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但 是 不 能 显示 表达 ， 数 值 计算 比较 耗 时 从 而 影响 配 准 速度 .寻找 二 元 Rayleigh 分 布 参数 更 优 的 
估计 和 二 元 Rayleigh 分 布 互 信息 的 近似 显示 表达 式 为 后 续 研 究 工作 . 
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Bivariate Rayleigh Distribution and its Mutual Information for SAR 
Image Registration 
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Abstract: Due to the speckle noise of Synthetic Aperture Radar (SAR) image which is determined 
by the SAR imaging mechanism, the interpolation artifacts of mutual information based on image 
registration become important. This paper combines the prior information that the homogenous regions 
of the single-look SAR amplitude image follow the Rayleigh distribution with the mutual information, 
and registers two single-look SAR amplitude images by the bivariate Rayleigh distribution and its 
mutual information. The experiments of synthetic images and real SAR images show the efficacy and 
advantage of this method when compared with the traditional mutual information method. 
Keywords: single-look amplitude; SAR; image registration; bivariate Rayleigh distribution; mutual 
information 
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